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Abstract 

We study the deviations from the exponential decay law, both in quantum field theory (QFT) 
and quantum mechanics (QM), for an unstable particle which can decay in (at least) two decay 
channels. After a review of general properties of non-exponential decay in QFT and QM, we 
evaluate in both cases the decay probability that the unstable particle decays in a given channel 
in the time interval between t and t + dt. An important quantity is the ratio of the probability 
of decay into the first and the second channel: this ratio is constant in the Breit-Wigner limit 
(in which the decay law is exponential) and equals the quantity ri/r2, where Fi and r2 are the 
respective tree-level decay widths. However, in the full treatment (both for QFT and QM) it is an 
oscillating function around the mean value Fi /F2 and the deviations from this mean value can be 
^Nj ■ sizable. Technically, we study the decay properties in QFT in the context of a superrenormalizable 

^ \ Lagrangian with scalar particles and in QM in the context of Lee Hamiltonians, which deliver 

CO . formally analogous expressions to the QFT case. 

(N 

^ 1 Introduction 

O; 

The study of unstable particles is a basic subject of both nuclear and particle physics. The fact that 
the so called survival probability, i.e. the probability p(t) that an unstable state S formed at the 
time t = did not decay at the instant t > 0, does not follow an exponential decay law for short 
times is now theoretically [TJ [H |31 HJ [51 15] and experimentally [71 |H1 IS] established in the framework of 
quantum mechanics. Moreover, oscillations on top of the exponential function can also occur in the 
short time regime and could eventually last long enough to be detected [TOl [11] . Also at very large 
, times deviations take place: a power law, and not an exponential law, is realized [1]. 

In particular, in Ref. [4] the deviation from the exponential decay law has been described in the 
nonrelativistic quantum mechanical context of Lee Hamiltonians [12[ I13j : this approach allows to 
mimic non-relativistic quantum field theory by including a continuum of states. In Refs. [141115] these 
deviations have been demonstrated in the framework of a genuinely (superrenormalizable) relativistic 
quantum field theoretical Lagrangian. Both approaches give rise to qualitative similar deviations from 
an exponential decay law in the short-time regime and also show an indiction of oscillations afterwards. 

In this work we study the decay of an unstable particle when more than one decay channel is present, 
both in the context of relativistic quantum field theory (QFT) and quantum mechanics (QM). This 
parallel allows also to show the formal similarities of these two approaches. The equations presented in 
this work are general and do not depend on the particular employed model(s). However, it is useful for 
a clear presentation to work with a well-defined theoretical setup: in QFT we employ, in line with Refs. 
|141ll5j. scalar fields, which are free from spin complications, in the framework of a superrenormalizable 
relativistic Lagrangian. In the case of QM we employ Lee Hamiltonians (LH): we generalize the work 
of Refs. ^ in order to describe the decay of an unstable state into two distinct channels. 
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When (at least) two decay channels are present, one has to determine (both in QFT and in QM) 
the probability of decay in each channel. Namely, the aim is to calculate (to our knowledge for the first 
time) the decay probability that the unstable particle S decays in the i-th channel between the time 
interval t and t + dt. Each partial decay probability deviates from the corresponding Breit-Wigner 
limit. In this context an important quantity under study in the two-channel case is the ratio between 
the decay probability in the first channel over the decay probability in the second channel. While 
this quantity is constant in the usual Breit-Wigner (BW) limit and is equal to the ratio T1/T2, where 
Fi and r2 are the tree-level decay widths in the first and the second channel respectively, this is not 
true in the general case. There are time intervals in which the ratio of probabilities is larger than 
ri/r2, meaning that the decay probability into the first channel is enhanced. Vice-versa, there are 
time intervals in which the ratio is smaller than ri/r2, implying that the decay probability into the 
second channel is enhanced. 

One of the main purposes of this paper is a general and theoretically based presentation of the 
equations related to the decay probability in each decay channel as function of time. For this reason 
we do not refer here to any particular physical system, but we keep the presentation as general as 
possible. Yet, in order to improve the understanding of the manuscript we shall present in the QFT 
case several plots of the introduced functions for a selected numerical example. Although the employed 
superrenormalizablc QFT model is clearly not the most general QFT treatment, the form an the 
structure of the equations are quite general and could be easily extended to other more complicated 
and realistic QFT approaches. 

The paper is organized as follows: In Sec. 2 we give an overview of the QFT formalism for the 
description of the decay of an unstable particle. We introduce a (superrenormalizablc) Lagrangian 
density and discuss the relevant properties of the survival probability p{t). In Sec. 3 we turn to the 
QM study of decays: to this end we introduce the Lee Hamiltonian. which shows a formally very 
similar behavior to the QFT case of Sec. 2. We review the most important equations in this approach 
and also present a relativistic generalization of it. In Sec. 4 we turn back to QFT and discuss in 
detail the formulae for the; two-channc;! case, with special attention to the decay probability in each 
channel. For sake of clarity we concentrate on the case in which two decay channels are available; 
the generalization to n decay channels is straightforward and is presented in Appendix B. In QFT 
the expressions of the decay probabilities in each channel cannot be analytically derived in the most 
general case; to circumvent this problem, a 'conjecture' about the form of the solution is put forward. 
It is then shown that the conjectured solution fulfills all the necessary requirements. A numerical 
study of these solutions is presented, in order to show differences from the usual exponential limit. In 
Sec. 5 we study the same problem of the two-channel decay in QM with the help of Lee Hamiltonians: 
in this case it is possible to derive the exact expressions for the decay probabilities in a given channel 
and we can test the validity of the conjecture presented in the QFT case: a numerical study shows a 
very good agreement of the exact results with the proposed solutions for both the nonrelativistic and 
relativistic treatments of the Lee Hamiltonians. Finally, in Sec. 6 we briefly present our conclusions 
and outlooks. 

2 General properties of decays in quantum field theory 

2.1 A simple Lagrangian to study the two-channel problem 

We introduce a relativistic (superrenormalizable) Lagrangian with the three scalar fields S, ipi and (p2- 

^ = \{d^Sr - ^M^S' + ^{d^^if - \ml^l + \{d,^2f - \ml<pl+grS^l +928^1 . (1) 

For definiteness, we assume that mi < 7712. The interaction terms giSLp\ and g2Sip2 induce the decay 
processes S (fifi and S — >■ (f2^2, see the pictorial representation in Fig. 1, panel (a). 
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It is useful to introduce the tree-level decay function as 



r'-'(x,TO,5) 



-{V2gf0{x -2m) , 



(2) 



where 0(x) is the step function and x plays the role of the 'running' mass of the unstable state S. The 
tree-level decay functions in the first and the second channel read respectively 



r^'\x,mi,gi) 



8ttx^ 



-(V2gife{x~2rni) 



''-{V2g2fe{x-2m2) 



(3) 



Similarly, we also introduce the total tree-level decay function as 

r's'(^) = r'<;^.^.(=^) + r*5U.W- (4) 

The partial tree-level decay rates -as calculated from the Lagrangian ([T])- are naively obtained by 
setting x = A/q into Eqs. ([3]). However, care is needed: as shown in the next subsection, the mass of 
the unstable particle S is modified by quantum fluctuationf0 (it is typically shifted to lower values): 



(5) 



For this reason the on-shell tree-level expressions Fi and F2 for the partial decay widths are not 
evaluated at Afp, but at the physical value M: 



Ti = = M) and F2 ^ Tf^,^^{x = M) 

Then, the total on-shell tree-level decay width reads 



F = F's'(a:-A/) = Fi+F2 



(6) 



(7) 



As a consequence, the tree-level (or Breit-Wigner) expression of the survival probability p{t) for the 
resonance S created at t = is the usual exponential function 



-Tt 



(8) 



and the corresponding tree-level expression of the mean lifetime is tbw = l/T- The explicit derivation 
of the exponential form is obtained through the Breit-Wigner limit in Sec. 2.3. 



2.2 The mass distribution ds{x) and the survival probability p{t) 

A crucial intermediate step toward the determination of the survival probability p{t) in the framework 
of quantum field theory is the evaluation of the propagator Gsip^) of the unstable resonance S. The 
quantity Gs{p^) is obtained at the 1-loop level by (re)summing the one-particle irreducible self-energy 
contribution (Fig. 1, panel (b)): 



Gs(/ 



- Ml + {^/2g^f^{p\ml) + (\/2g2)'S(p', m^) + 



le 



(9) 



where T,(jP,m^) is the amplitude of the standard loop diagram in which two virtual particles of the 
(fi-type circulate: 



1 



{2nY 



q) — + te 



(I — g) —171^ + is 



(10) 



^Alternatively, one could introduce a counterterm and work with Mq as the physical mass, see details in Appendix 1. 
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(a) Decay (bj Seif-energy 

Figure 1: Pictorial representation of (a) the amplitude of the tree- level decay process S" — >■ ^pitpi with 
i — 1,2 and (b) the corresponding self-energy diagram. The imaginary part of (b) is related to (a), 
see Eq. (IT3l) and the discussion in the text. The quantity (p{q) is the vertex function, which is chosen 
to be 4i{q) = 9{A? — q^) for the numerical evaluations. 



It is also useful for future purposes to introduce the function n(p^) as 

n(p2) = (V2gif^{p\ml) + (V2g2)'S(/,m2) 
in such a way that the propagator is given by the simple expression: 



- le 



(11) 



(12) 



Due to the optical theorem the imaginary part of the function Ii{p^ — x^) is related to the tree-level 
decay function in Eq. 

Imn(a;2) = a;r*/(x) . (13) 

Although the imaginary part is convergent, this is not the case for the real part. For this reason the 
loop function I](p^,to^) must be regularized: in this work we use a sharp, three-dimensional cutoff A 
for the explicit numerical examples. However, different choices of the regularization procedures do not 
change the qualitative features of the discussion [T31 [HI [H] . The explicit expression for the function 
S(p^,to^) and more details on the consequences of the optical theorem are given in Appendix A. 
The spectral function ds{x) of the scalar field S is defined as the imaginary part of the propagator: 



ds{x 



V) = — 
TT 



Explicitly: 



ds{x) 



2x 

— lim 

TT 6^0 (a.2 



Im [U{x^)] 



M§ + Ren(a;2))- -f- (Imn(a;2) + e)" 



The renormalized mass M of the particle S is defined as the zero of the real part of Gs{p^) ^ 

- Ml -f Ren(M2) = . 
Usually, Ren(Af2) > o, thus the renormalized mass M is smaller than Mq. 



(14) 



(15) 



(16) 
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Figure 2: The spectral function ds{x) from Eq. (fT4|) and the nonrelativistic BW limit ds.Bwl^;), see 
Eq. (pOl derived in Sec. 2.3, are plotted for the numerical values of Eq. ([TSl) . In this example the 
function ds{x) differs only slightly from the BW limit, but this small deviation is enough to generate 
sizable deviations in the early temporal evolution, see later on. 



The quantity ds{x)dx represents the probability that in the rest frame of S the state S has a mass 
between x and x + dx. For this reason it is also denoted as the mass distribution of the unstable state 
S. It is correctly normalized for each 51,52, 

/•oo 

/ dxdsix) = 1 , (17) 
Jo 

and, in the limit 51 — > and 52 0, the expected spectral function ds{x) = S{x — AIq) is obtained, 
see details in Refs. [TBI [T7] . 

In Fig. 2 we present the spectral function ds{x) for a selected numerical example, where we chose 
as energy unit the quantity [2m2]: 

Mo = ^ [2m2] , 31 = .92 - 3 [2to2] ' '"i = ^ [2™2] : "I2 = ^ [2m2] , A = 2 [2m2] . (18) 

As a result one obtains from Eq. (fT6)) that M — 1.40 [2to2] (which is just slightly smaller than 
Mq: the quantum fluctuations are not large in this example). The tree- level decay widths read: 
Fi = 0.13 [2m2], F2 = 0.089 [2m2] and F = 0.22 [2m2] . As a consequence, the tree-level lifetime reads 
■''BW — 4.63 [2m2] . We shall use these numerical values in all the upcoming plots of the paper, 
because they allow for a clear visualization of the underlying equations. It should be however stressed 
that the qualitative behavior is not dependent on the precise numerical choice. 

It is now possible to determine the probability amplitude a(t), and therefore the survival probability 
p{t). To this end, let us consider the set {\x)} of eigenstates of the full Hamiltonian H of the system 
defined by the Lagrangian of Eq. ([T]): H \ x) — x\x) . The continuous variable x describes two-particle 
states and therefore is such that x > 2mi. The survival amplitude a{t) reads (in the Heisenberg picture) 

a{t) ^ {S\e''^'\S) = {S\x) e-'""* {x\S) = |(5|x)|'e-"* , (19) 

where the integration over the continuos variable x is understood. Taking into account that |(S'|x)|^ = 
ds{x) (i.e., as shown above, ds{x)dx is the probability that the unstable state S has an energy between 
X and X + dx) and making the integration explicit, we get: 
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ait) - / dxd5(a;)e~"* , p{t) = \a{t)\^ . (20) 



The integrand in Eq. ([20|) is nonzero only between (2mi, oo) in virtue of the step function in Eq. (fT3|) . 
However, it is from a technical point of view helpful to write the integration range as (—00, 00) so that 
the quantity o(t) is expressed as the Fourier transform of the spectral function ds{x). It is in this way 
that the usual exponential law can be derived in the so-called Breit-Wigner limit, see Sec. 2.3. 

The condition p{0) = 1 is fulfilled in virtue of the normalization of ds{x). This property is, in turn, 
a consequence of the resummation and of the validity of the so-called Kallen-Lehman representation 
[IT] . In Fig. 3 the function p{t) and the corresponding exponential limit pBw{t) — e~^* are plotted: it 
is evident that sizable deviations from the 'full' result p{t) take plac(3. This fact is remarkable if one 
considers that the deviations of the spectral function ds{x) from the Breit-Wigner limit ds-Bw{x) are 
very small, see Fig. 2. Nevertheless, these small deviations are enough to generate large changes in the 
temporal evolution of the system: for a better visualization of this point, the difference p{t) — pbw(0 
has been also plotted in Fig. 3. 

The spectral function ds{x) is naturally bound from below because of the presence of the low- 
energy threshold 2mi. Its behavior for large energies x can be easily inferred by using the optical 
theorem in Eq. p^ : ds{x) ^ x~^ for large x. Then, the presence of cutoff A implies that ds{x) — 
for X > 2-\/A2 + mf ~ 2A. It thus follows that the mean energy 

dxxdsix) (21) 

) 

is a finite number. As shown in Ref. [I] the existence of an energy threshold and the finiteness of 
(x) suffice to guarantee short-time (and also and long-time) deviations from the exponential decay 
law. Quite remarkably, in the present superrenormalizable model of Eq. ^ the finiteness of (x) is 
guaranteed also in the limit A — > 00, being xds{x) ^ for large x. On the contrary, all higher 
momenta (1") = J^^dxx"ds{x) with n = 2, 3, ...diverge for A — 00. As shown in Ref. [TH], the 
finiteness of (x) also assures the differentiability of the amplitude a{t) (and thus also of p{t)), which 
we will use in the following considerations and also in Sees. 2.4 and 2.5 when discussing in more detail 
the short-time behavior and the quantum-Zeno effect. 

We now introduce the decay probability density h(t): the quantity h(t)dt represents the probability 
that the unstable particle S decays in the time interval (t.t + dt). Clearly, the integrated quantity w(t) 
defined as ^ 

w{t) = [ A\xh{u) = l-p{t) (22) 



Jo 

represents the probability that the unstable particle decays between and t. Deriving the latter 
equation we obtain 

Mt) = -f = -p'W , (23) 

i.e. h(t) is simply the negative derivative of the survival probability p{t). 

In terms of state vectors the probability w{t) can be symbolically expressed as 



,{t) = f dnh{u) = ^ |(¥'^.k^^,-k \e-'"'\ S)\ 

•^0 t=l,2 k 



(24) 



where the state |((5i,k</'i,-k) represents the state of two particles of the type tpi, which carry opposite 
three-momentum k and — k, respectively. In fact, the object | (</'i,k<^i,-k |e~*^*| is the probability 
that the particle S decays into two particles of the type ipi, the first with momentum k and the second 



^Here 'full' result means that the resummed 1-loop approximation has been performed and that the full functional 
form of the spectral function ds{x) is kept in the evaluation of the temporal behavior. 
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Figure 3: The survival probability p{t), the BW limit pBw(i) = e~^* and the difference p{t) — PBw{t) 
are plotted for the numerical values of Eq. ([TSl) . The dot on the horitonzal axis corresponds to the 
value of tm, see Eq. (^0)1 and the discussion in Sec. 2.4. 



with momentum — k. The overall decay probability w{t) results as the sum of all these (independent) 
probabilities. It is important to stress that in Eq. ([Ml) we neglect the decays into n-particle states 
with n = 4, 6, 8, .., which are suppressed by kinematical factors and anyhow neglected in the self-energy 
contributions. Note also that the summation means the sum over all the momenta k = 27rn/L, 
where n = (ni, 712,^3) and rii = 0,±1,±2,...; the system is thus confined into a box with volume 
V — L^. The limit L — > cxo is straightforward but unnecessary for the present discussion (see Sec. 3.1 
for the explicit evaluation of it). 

Upon deriving Eq. ([M|) we can also express h{t) in terms of state vectors: 

hit) = jAY.T.\(^^M^^^-^\e-'"'\S)n . (25) 

yi=l,2 k J 

It is indeed a remarkable fact that, in order to evaluate h{t) we do not need to evaluate explicitly the 
matrix elements (</'i,k<Pi,-k |e~*^*| S*), but we simply have to perform the derivative of Eq. (I23p . 

In Fig. 4 the density of decay probability h{t) and its corresponding BW limit /ibw(0 — — Pbw(^) — 
Fe"'"* are plotted. Here the difference between the 'full' result and the BW limit are even more evident, 
due to the fact that h{t — > 0) vanishes and h-Qwit — >■ 0) = F clearly does not, see Sec. 2.4 for a detailed 
description of this point. Moreover, one can clearly see the presence of oscillations of the function h(t), 
in agreement with the discussion of Ref . [TT] . 

2.3 The Breit-Wigner limit 

It is instructive to discuss the Brcit-Wigner limit, in which the survival probability p{t) turns out to 
be exponential. We start by approximating ds{x) as follows: 

M^) = - '-^^ 2M MT^,HM)eix-2m,) 

^ {x^- M§ + ReUix^)f + {xT's\x)f ^ (x^ - M^f + {MTf{M)Y 
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Figure 4: The decay probability density h{t) and the corresponding BW hmit /iBw(i) = Fe^'"* are 
plotted for the numerical choice of Eq. ([TSl) . For short times the two functions are qualitatively (and 
here also quantitatively) very different, while the former vanishes and the latter does not. Moreover, 
oscillations on top of the exponential are visible. Then, for large times the function h(t) approaches 
the BW hmit /iBw(i)- 



By having included the step function 9{x — 2mi) we have kept track of the threshold. A subtle point 
of Eq. (|26p is the normalization, which is spoiled by the employed approximation. This property is 
however crucial in order to have p{0) = 1. It is thus necessary to rescale the function in such a way 
that it is normalized to 1. Finally, the spectral function in the relativistic BW approximation reads: 

-^BWrel 



ds-BWrel{x) 



(a;2 - Ai"2)2 + iMVy 



where 



BWrcl 



dx- 



e{x-2mi) 



-\ -1 



-e{x - 2mi) 



dx- 



1 



(x2 - Apy + {Mvy 

A further simplification, leading to the usual, nonrelativistic Breit-Wigner form, is needed: 

A^BWrci6'(a; - 2toi) iVBWroi6'(a; - 2mi) 



(27) 



(28) 



ds- 



BWrcl 



(^) 



{x^ - Apy + {Airy 

Nbwic\0{x ~ 2mi) 



{x ~ Aiy {x + Aiy + {Airy 

A^BWrol 



■ (29) 

{x ~ Aiy {2Aiy + {Mry {x - Aiy {2Aiy + {Airy 

Note that in the last equation we have also dropped the requirement that x > 2mi, thus x can now 
vary from — oo to +oo. This is of course unphysical for the present Hamiltonian because negative values 
of X, i.e. negative values of the running mass of the unstable state S, do not make sense. Nevertheless, 
the function is peaked at M > and the extension to negative x affects the function for x > 2toi only 
negligibljl^. As a result, the spectral function in the nonrelativistic BW approximation reads 



•^This is true for the nonrelativistic approximation, but is not true for the relativistic one: in fact, the latter is peaked 
also for X = —M. 
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where -/Vbw = r/27r is the normahzation constant such that dxds-Bwix) = 1. 
The aniphtude a{t) as calculated from Eqs. (|20p and ([501 takes the form 

aBw(0 = e-"/2e-^^^* {t > 0) (31) 

which implies the exponential decay law 

PBwW = laewWl" =e-" . (32) 

Notice that this result has been obtained by extending the range of integration in Eq. ((20)) to (— oo, oo). 
Obviously, the Breit-Wigner distribution is never exact and cannot be fulfilled by any physical system. 
It is in this sense remarkable that deviations from the exponential decay law have been discussed 
only relatively recently. This has surely to do with the fact that the BW distribution, although never 
exact, works extremely well for the majority of physical systems in a sufficiently large energy interval. 
Conversely, as function of time the deviations from the exponential law take place very early and very 
late, and are thus difficult to measure. 

The decay probability density h{t) defined in Eq. (|23|) shows also an exponential behavior in the 
BW hmit: 

h{t) -> hnwit) = -PBw(i) = Te-r* . (33) 

The exponential functions PBw{t) and h^wit) have been plotted and compared with the 'full' results 
p{t) and h{t) in Fig. 3 and 4, respectively. 



2.4 Vanishing of p'{0) and the definition of the 'Zeno-time' 

Let us consider the first derivative of the survival probability p{t) for t — 0, p'(0). To this end we 
study the quantity 

POO 

a'{t) = / dx(-i.T)ds(a;)e~"* . (34) 



In virtue of the cutoff (independently on its value) the integral dxxds{x)e is well defined and 
does not contain any divergences. In particular 

/oo 
dxxdsix) = -i{x) (35) 
-00 

is a well defined imaginary number, being (x) a finite real number, as discussed in Sec. 2.2. 
The function p'{t) can be expressed as 

p'{t) = a'*{t)a{t) + a*{t)a'{t) . (36) 

Being a(0) = 1 we obtain 

p'(0) = i (a;) - j {x) = 0. (37) 

The first derivative of the survival probability vanishes for t = 0. This means that for small times a 
deviation from the exponential law, for which one obviously has limj_jo+ P'swi^) ~ (^^^ -^Q- (EH))i 
is realized. 

The natural question is how long does the deviation form the exponential law lasts. In the literature 
a Taylor expansion has been often performed, 

p{t) - 1 - 4 + - : (38) 

where the coefficient 

- - 
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has been defined as the 'Zeno time', which quantifies the time interval of non-exponential behavior. 
However, as already explained in Ref. |14| . the second order derivative does not even need to exist. 
As a trivial illustrative example one can consider the simple and well defined function p[t) = 1 — t^^'^ , 
for which p'{0) = but p(">i)(0) = oo. Moreover, even when the second derivative exists -as in all 
physical cases, see below- this does not necessarily mean that the parameter tz in Eq. (pS)) represents 
a correct estimate of the time interval in which a deviation from the exponential law takes place. 

In order to overcome this problem, in Ref. |14| a definition has been introduced which does not 
depend on the Taylor expansion of the function p(t) . The time tm is defined as the instant of time at 
which the deviation of the function p{t) from the exponential behavior e~'"* is maximal: 

max {p{t) - e"^*) t ^ tm ■ (40) 

Clearly 

= (41) 

in all practical cases tm corresponds to the first zero of the derivative of the function p{t) — e^'"*. The 
function p{t) — e~^* has been plotted in Fig. 3 for the numerical values of Eq. (1181) . The maximum 
is realized for tm = 1-18 [2m2] , implying the ratio tm/tbw = 0.26. In the present model, the non- 
exponential behavior elapses for ~ 2tm, meaning that in this example is about 50% of the tree-level 
lifetime of the unstable particle. 
Further comments are in order: 

(i) In the present work we consider field theoretical examples with some finite value of the cutoff. It 
is indeed a natural requirement that at some energy x > A the loop function I](a;^,m^) is suppressed. 
The actual value of the constant A depends on the theory under study. For instance, working with 
hadronic theories implies A ~ 1-2 GeV [Hj. On the contrary, for elementary Standard Model processes, 
the cutoff A is of the order of Mpianck (alternatives are A ~ Mqut ~ 10^^ GeV and, in some recent 
scenarios with large extradimension(s), even smaller with A > 1 TeV). 

(ii) One can also regard the theory of Eq. ([T]) as a (super)renormalizable toy model with a very 
large value of the cutoff, such as A ~ Mpianck 10^^ GeV. In this case the first derivative vanishes, 
p'(0) = 0, and all the other derivatives, including the second derivative p"(0), are large but finite: for 
instance, the quantity p"(0) is proportional to In A. This is easily understandable, being the asymptotic 
behavior for large x given by ds{x) ~ , thus implying that (a;^) = ^_^x^ds{x)Ayi ~ In A. In the 
present superrenormalizable case, the functions p(t) for A ^ 1-2 GeV and for A ^ Afpianck differ only 
slightly: the deviations from the exponential law last for a similar amount of time and the value of the 
ratio tm/tbw shows a very weak dependence on the cutoff [TJ. This example shows that, in the case 
of large A (as for instance, A ^ Mpianck), although the Zeno time tz ~ 1 / Vln A turns out to be very 
small (as it was obtained in Refs. [20l[2T] in the framework of perturbation theory), tz offers only an 
underestimation of the actual non-exponential time interval, which is sizable also for large values of 
the cutoff. 

(iii) In the nonrelativistic Breit-Wigner limit of Eq. (l30l) the first derivative p'(0) is -strictly 
speaking- not defined. However, by using the fact that the amplitude reads pbw(^ > 0) = e~'"*, one 
obviously has lim(^o+ Pbw(^) = other words, the function PBw(i) can be easily continued to 
the point t = 0. 

(iv) The present results are based on the (resummed) 1-loop approximation. The inclusion of 
higher order, such as the sunset diagram in which a particle S is exchanged by the two particles if 
circulating in the loop of Fig. l.b, is a task for the future. It is expected that higher order terms do not 
change substantially the results: in fact, higher order diagrams are suppressed by the multiplication 
of vertex functions. (Moreover, in the framework of hadronic theories, higher order are suppressed in 
the so-called large- A^c approximation [3^. The behavior of a theory similar to that of Eq. ([T]) in the 
large- iVc limit is discussed in Ref. [53]). 

(v) The resummation of the loop diagram in Fig. l.b, expressed in its unregularized form in 
Eq. (|10|) and in its regularized form in Eq. ()144|) reported in the Appendix A, has been a crucial 
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step toward the presented results, also in connection with the correct normalization of p{t) and the 
vanishing of p'(0) (necessary to show the non-exponential behavior). If, instead, one would consider 
only the perturbative expressions of the propagator at the 1-loop (gf) level, one obtains 

GT'ip') = , L , + ( , L , - V \{V2g,mp^ml) + (^2.92)^^(^2, m^)! + ... , (42) 



thus the spectral function in this case would take the form 



d, (.) = -ImG| ip ^.)^5{x~ Mo) + — _ + V (x^ - Mg)^ + - ^^^^ 

The normalization of c?^'''^'(a;) is not anymore explicitly fulfilled for each value of g and an unphysical 
singularity for the bare value of the mass x = Mq is present. (This infinity is removed when the 
resummation is performed.) It would be cumbersome to derive analytical results with this perturbative 
expressions; moreover, a numerical study could hardly be performed, thus showing the necessity of the 
resummation before the temporal analysis is investigated. 

(vi) The issue of the preparation of the state \S) at the time i = is well defined when the state is 
long lived, but more controversial for short living particles. A related difficulty is the (in) dependence 
on the representation choice [21] , see also the corresponding discussions in Ref. [13] . Here it sufiices to 
say that the change of the representation of the fields is equivalent to the change of the initial state. 
Nevertheless, also in this case deviations from the exponential law take place, but the mass distribution 
would have a different form. 

(vii) As a last step we discuss the validity of the employed model of Eq. ([T]) upon variation of the 
parameters. The choice of the numerical values in Eq. ()18p has been done in order to obtain figures 
which allow for a clear presentation of the underlying equations. However, the described patterns 
take place also when the parameters are varied: in particular, when the coupling constants gi and 
g2 are decreased, the time interval in which sizable variations from the exponential decay law take 
place becomes smaller. In the opposite direction, when the coupling constants gi and 52 become 
larger, this time interval grows. However, care is needed when increasing the values of the coupling 
constants: beyond a certain threshold, the validity of the Kallen-Lehman representation, and thus 
of the normalization condition of Eq. ([T7)) . does not any longer hold, see the detailed discussion in 
Ref. [TTj. Note also that the QFT model of Eq. ([ij does not show, for any value of the parameters, 
repopulation effect. In fact, the function h{t) is always positive in this model. 

2.5 The quantum Zeno effect 

It is interesting and amusing to briefly review the so-called quantum Zeno effect. To this end we deflne 
the function j{t) as in Ref. [3]: 

^(i) = _i!^op(t)-e-^(*)* . (44) 

In Fig. 5 the function ^{t) is plotted and is compared to the tree-level decay width F. It is visible that 
for short times "f(t) is considerably smaller than F, see the discussion below. 

Let us now suppose to perform a single measurement of the system at the time t > 0. The 
probability that the unstable state S did not decay at the time t reads 

= p(t) = e-'^(^)^ ~ e-^^ (45) 

where in the last passage we assumed that the time t is large enough in order that 7(t) ~ F. 

Instead, let us suppose to perform N measurements at the time intervals t*, such that r = Nt^,. 
The probability that at the time r (after the N measurements) the state did not decay is given by 



PN-t. = (e-^(*-)**) = e-^(**)^** = e-^(*-)^ . (46) 
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Figure 5: The function ^{t) = — [Inp(i)] /t is plotted. For comparison, the constant BW limit 7bw(0 = 
— [InpBw(^)] /t = T \s also shown. For short times ^{t) < F, rendering the quantum Zeno effect 
possible. As usual, the numerical values of Eq. ([T8| have been used. 



When belongs to the 'non-exponential' time-interval, < < tm, one has 7(t*) < F. Indeed, 
considering that the function ^(t) vanishes for t = 0, 



holds. Its meaning is clear but astonishing: the subsequent measurements of the unstable particle 
S increased the probability that it did not decay. (Formally, one performs a 'measurement' of the 
unperturbed Hamiltonian Hq, in which the coupling constants gi are set to zero). This property has 
been indeed verified experimentally by using ions trapped in a potential barrier [9]. 

For — > 0, i.e. for iV — cx) (a situation which corresponds to a 'continuos' measurement) one 
obtains that pN-t, — ^ 1- The unstable particle does not decay at all! This is the famous quantum Zeno 
paradox. 

As Fig. 5 shows, there are also time intervals in which 'y{t) > F. This implies that for satisfying 
this relation a faster decay rate is obtained when the state of the unstable particle is observed at 
intervals of . This situation corresponds to the so called quantum Anti-Zeno effect |S1 H] ■ 

We conclude this brief recall of the Zeno effect by noticing that it does not take place in the 
Breit-Wigner limit. In fact, in this limit one has j{t) — F for each thus meaning that pN-t, = 
Pt, independently on the value of and t. Making a single measurements or many intermediate 
measurements does not lead to any change in the survival probability. 




(47) 



one can make the quantity 7(i*) as small as desired by choosing ^tm- 
In general, for < < 2tm, the inequality 

PN-t, ^ e-''^''^^ > Pr c^e-r^ 



(48) 
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3 General properties of decays in QM: Lee Hamiltonian and 
link to QFT 

3.1 QM with Lee Hamiltonian(s) 

We now study a similar system in the framework of quantum mechanics. In order to obtain irreversible 
dynamics in a QM approach, one needs to couple the unstable state \S) to a continuum of states. The 
formalism described here relies on the so-called Lee Hamiltonians [12 (later on denoted also as QM-LH 
approach) and is a generalization to the two-channel case of the work in Ref. [4 . In this section we 
recall general properties of the LH approach with the special aim of showing the similarities with the 
genuine QFT approach presented in the previous section. In this context we also refer to Ref. jl3j . 
where a solvable Lee Hamiltonian has been studied in great detail and the relation to QFT has been 
described. 

For transparency, we first present the Hamiltonian using a discrete sum and later on perform the 
transition to the continuous results. The Hamiltonian under consideration is given by 

H = Ho + H^nt , = + H^, , (49) 

where 

ffo-Mo|5)(5|+^c^i(k)|l,k)(l,k|+^t^2(k)|2,k)(2,k| , (50) 

k k 



ij£)=^gi^(|l,k)(5|-|-|5)(l,k|) , (51) 
iji^)=5]52^(|2,k)(5| + |5)(2,k|) . (52) 

k 

In order to establish a contact with the study of Sec. 2, the variable k stays for k ~ 27rn/L and 
V — Li^. The analogy is clear: the state 15) corresponds to the particle S of Eq. ^ and the states 
|l,k), |2,k) are analogous to the two-particles states Itpi.kVi.-k) and |</32,k</'2,~k) of Sec. 2. The mixing 
terms in Eqs. ((5T|) - ([52)) correspond to the interaction terms giSip\ and g2S(p2 of Eq. JT]), respectively. 
The transition to the continuous limit is obtained for 



T.'^^J^^\^^^)^r^\^'^) = (53) 
The time-evolution operator is given U{t) = e~'^* can be expressed in terms of a Fourier transform: 
U{t) = e-^' = - / dE rr^. e-^"' = - / dEG{E)e-^^' , (54) 
where the operator 

has been introduced. One can rewrite G[E) as follow: 

1 ^ 1 1 _ 1 ^ / 1 \" 

^ E-H + ie " E-Ho + iel- ^ ^ E-H^ + ie \^'"' E - Ho + ie ) 

The propagator Gs{E) of the unstable state \S) is defined as 

Gs{E)^{S\G{E)\S) (57) 
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and can be evaluated using Eq. ([56| : 



E — Mq + IS ^-^ \ E — Hq + le I 

n—O 



Explicitly, the cases n = 1 and n = 2 read: 



n = l^^-^(5|F„.|^)=0, (59) 

The recursive quantity found for each n even is Il{E): 

U{E) - - (5| ^ ^ . 1 5) = (61) 

- (5| (|1, k) (1, k| + |2, k) (2, k|) ^ (|1, q) (1, q| + |2, q) (2, q|) |5) , (62) 

where the sums over k and q are understood. Using Eqs. ([5T|) - ([52l) and the orthonorniality of the 
basis we get: 

k k 

In the continuous limit the basic quantity n(_E) becomes: 

d'k gfffjk) r d^k ff|/|(k) 

{2n)^ E - uJi{k) + ie J (27r)3 ^ - W2(k) + ie 

The propagator Gs{E) of the state |5) reads then 

Gs(E) ^ {S\G(E)\S) = — ^— . (65) 

The similarities between the present QM-LH approach and the QFT approach described in the previous 
section are clear (see for instance the form of the propagator in Eq. the main difference being 

the emergence of nonrelativistic propagator in the present case. The quantity Tl{E) represents the 
analogous of the loop contributions. 

We can now calculate the survival amplitude a{t) = {S\ U{t) \ S): 

a{t) = {S\ U{t) \S) = {S\- dE———-e~^^' 1^) = / dEGs{E)e~^^K (66) 

27rJ„oo E-H + ie ^t^J-oo 

Using the dispersion relation 

T^J-oo x-E-ie 
the survival amplitude a{t) can be expressed as 

/oo 
dxd5(a;)e-"* (68) 
-OO 

where in this context the identification 

ds{x)^-\lm[Gs{x)]\ (69) 

TT 

has been made. Again, the expression and its QFT-counterpart of Eq. (IT?)) are formally equivalent 
(a part from unimportant normalizations). Obviously, all the discussion about the deviations from the 
exponential law, the Breit-Wigner limit, the vanishing of p'(0) and the quantum Zcno effect can be 
easily repeated in the present QM-LH context. 
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3.2 Relativistic generalization of the LH and link to QFT 



In order to obtain expressions which are equal to those of QFT of Sec. 2, one has to consider a 
slight modification of the Lee model. The issue is that in QFT one has particles with positive energy 
propagating forward in time and particles with negative energy propagating backwards in time (i.e. 
antiparticles). In our QFT model of Eq. ([Ij particles and antiparticles coincide, but it is anyhow 
the combination of both contributions which give rise to the usual relativistic propagator of the form 
— TO^ + ie). 

The QFT expressions can be obtained in a LH approach by modifying the time evolution operator 
([54]) in the following way: 



U{t) = ^ dE ( ^ I . + ^ I . ) e-'"' (70) 

^' 2ttJ_^ \E-H + ie E + H-ieJ ^ ^ 

+ 00 njTi ■ p + OO 



■ 2E 



= — dE 



-iEt 



2n E^ - 

where the operator G{E) reads now 



le 



— / dE2EGiE)e-'^\ (71) 



The integral has been done according to the famous Feynman prescription. If, as a simple example, 
one considers a free particle \S) for which Hq = Mq \S) {S\, one finds the usual propagator form (in 
the rest frame of S) GsiE) = {S\ G{E) \S)^ll{E'^~ + ie). 
We consider a model in which we can write as 

= + HI^^, Hf^^ H^^J^ + h\^^ (73) 

where 

Ho'=Mo l^)(5|+E^iW|l'k>(l.k|+E^2(k)|2,k)(2,k| , (74) 



ijW^ = X]5i^(|l,k)(5| + |5)(l,k|) , (75) 
^LT = E52^(|2,k)(5| + |5)(2,k|) . (76) 



j(2)2 _ „ /2(k) 

k 

We then follow the same steps of Sec. 3.1 and introduce the propagator as 



GsiE) = (^1 GiE) 15) = ^..^.^^(i.)^,^ , (77) 

where n(i?) is calculated in the same way as in Sec. 3.1 and -as expected- turns out to be the 
relativistic generalization of Eq. ([M)) : 

f d^k g'ifK^) f d^k gi/|(k) 9 , , ■? , . 

J (27r)3 — w^(k) + z£ J (27r)3 ii^ — a;2(k) + 

Notice that the propagator Gs{E) = {S\ G{E) \S) takes a form which is formally equivalent to Eq. ^ 
obtained in QFT. 

In particular, when we make the identifications 

c^?(k) = 4(k2 + m\) , c^2(k) = 4(k2 + ml) , (79) 
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and 



/i(k) = V2- 



0(k) 



we get indeed in the present Lee model the same mathematical expression for the propagator of the 
QFT approach of Sec. 2. Note, the additional factor \/2 is here necessary to make contact with the 
identical particles in QFT and the cutoff function (/)(q) is described in detail in Appendix A and takes 
the form (i>{c\) = ^(A^ — q^) for the numerical evaluations presented in Figs. 2-10. Loosely speaking, the 
QFT case corresponds to a particular choice of the functions Wi(k) and /i(k) of the Lee Hamiltonian 
studied here. This is an interesting property which we will use in Sec. 5. 

By making use of dispersion relations we arrive at the usual expression for the survival amplitude 
a{t): 



, /2(k) = V2- 



m 

ml) 



1/4 



(80) 



a{t) — I dx(is(a;)e 



where 



ds{x) 



2x 

TT 



lmiIm[(S'|Gs(x)|5)] 



(81) 



(82) 



is identical to Eq. (|14|) . This shows once more the tight connection of the present formalism with that 
of Sec. 2. 



4 Two-channel case: quantum field theory 

4.1 The decay probabilities Wi{t) and W2{t) and the corresponding densities 
hiit) and h2{t): the formal expressions 

In this section we come back to the quantum field theoretical model of Eq. ([T]). The central quantities 
under investigation are the decay probability densities hi (t) and /12 (t) and the corresponding integrated 
decay probabilities wi(t) and W2it). The object hi(t)dt represents the probability that the particle S 
decays in the first channel {S — > fi(pi) in the time interval between t and t + dt; as a consequence, the 
integral wi{t) — du/ii(u) represents the probability that the particle S decays into the first channel 
in the time interval {0,t). The quantity wiit) can be formally expressed as the sum of the probabilities 
of decays into ififi for different values of the three- momentum k: 

wi{t)= / du/ii(«) = V|((^i,u(^i,_k|e-'^*|5)|' . (83) 
Jo k 

Similarly, h2(t)dt represents the probability that the particle S decays in the second channel {S — >■ 
ip2f2) in the time interval between t and t + dt and the integral 

W2{t) = f duh2{u) = V |((^2,k<^2,-k \e-'"'\S)\^ (84) 
Jo k 

represents the probability that the particle S decays into the second channel in the time interval (0, t). 
The equality w{t) — wi{t) + W2{t) holds by construction (see Eq. (p4)) ) and is intuitively clear: the 
total decay probability between (0, t) is the sum of the decay probabilities in each channel in the same 
time interval. 
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The symbolic expressions for hi{t) and ft.2(t) in terms of state vectors are obtained by deriving Eqs. 
and (jM)) : 

hi{t) = i{y:\(vi,uVi,-^\e-'"'\S)r] , (85) 



dt 



u(p2,-k|e-'^*|5)|'j . (86) 

In this section the decay probabihty densities hi{t) and /i2(i) are not calculated directly, because 
the evaluation of the matrix elements ('ySi.kV'i.-k |e~*^*| 5) and the subsequent partial sum is not an 
easy task in QFT. We will follow a different way by starting from the decay probability density h{t), 
which by construction is the sum of hi{t) and h2{t). It is then possible to postulate the expressions 
for the latter quantities and then a posteriori verify that they fulfill all the required properties. We 
will check the numerical validity of our solutions with the help of the QM-LH approach (both in the 
nonrelativistic and the relativistic frameworks) in Sec. 5. 

An important object of interest in this section is the ratio of decay probability densities R{t) defined 

as 

It represents the ratio of the probabilities that the particle S decays in the first channel and in the 
second channel in the time interval between t and t + dt. Its study is interesting because, as we shall 
see, R{t) shows a peculiar behavior as function of time. 

For a better understanding of the presented functions several plots for the numerical values of Eq. 
(lis)) are presented. 



4.2 The partial distributions Sg\x) and d^g\x) and amplitudes ai{t) and 
The distribution ds{x) can be written as the sum of the two functions d^g\x) and d^g\x): 



with 



dsix)=d^p{x)+S^\x) , 

2x 2gllm [J:{x^,ml)] +e 

da '(X) = — hm 5 , (89) 

TT e->o (2-2 _ ^2 ^ Ren(a;2))^ + (Imn(2;2) + ef 

4^)(x) = ^ lim 2gllm[nx\rni)]+e 

^ (x2 - M2 + Ren(a;2))' + (Imn(a;2) +£)' 



The integrals 



/>oo />oo 

/ dx4^^(a;) and / dxd^s\x) (91) 
Jo Jo 

represent the decay fractions, i.e. the probabilities that the particle S decays in channel 1 or in channel 
2 in the time interval (0,oo), respectively |25j . Obviously, their sum is unity in virtue of Eq. (jl7p . 
Note, in the BW limit the integral (hid^g\x) represents the ratio Fi/F, see Sec. 4.3.5 for details. 

The functions d^g\x) and d^g\x) have been plotted in Fig. 6 for the numerical values of Eq. p^ . 

We now introduce the Fourier transforms of d^g\x) and d^g \x): 

/>oo 

ai{t) = / dxd^g\x)e-"'' , (92) 



17 



a2{t) = / dx4^^(a;)e-*^* 



The survival amplitude a{t) can be rewritten as 



a{t) = / dxds{x)e 



dx4^^(a;)e-"* + / dx4^(a;)e'"* = ai{t) + 03 (t) 



Then, the survival probability p{t) can be expressed as: 



p{t) = a*{t)a{t) = (alit) + al{t)) {a^{t) + 02^) 

= a\{t)ai{t) + a\{t)a2{t) + ai{t)a*2{t) + a2{t)a2{t) 

= Ai{t)+2ArmAt)+A2{t) , 



where the three real quantities 



A,{t)^\a,{t)r 



A2{t) = \a2{t)\' = 



dx^-* (x)e 



Ajnix (i) 



al{t)a2{t) + aiit)a*2{t) 



Re[ai{t)a2it)] 



have been introduced. 

Therefore, the decay probability density h{t) — ~p'{t) can be rewritten as: 

h{t)^-A[{t)~2A[^,,{t)^A'2{t) . 



(93) 



(94) 



(95) 

(96) 
(97) 
(98) 

(99) 



This form is the starting point for the postulate about hi(t) and h2{t) as presented in the next 
subsection. 



4.3 The decay probabilities Wi{t) and 1^2 (t) and the densities hi{t) and h2{t) 
as functions of ai{t) and a2{t): a conjectured solution and its properties 

4.3.1 The postulate about hi{t) and 'Wi{t) 

The key functions in the present study are the densities of decay probability per each channel, denoted 
as hi{t) and h2{t), which have been introduced in the formal expressions (f85|) and (f86|) . 
We postulate that the functions hi(t) and /i2(i) can be calculated as: 

h,{t) ^ ~A[{t) ~ A'^,,{t) , (100) 

h2{t) = -A'2{t) - A'„,,,{t) . (101) 

The direct calculation of Eqs. ([85]) and (|86l) is a complicated task in QFT, see also next Section where 
the explicit evaluation in terms of Lee Hamiltonian(s) is performed. Here we argue that it is possible 
to express the final result by using the introduced partial amplitudes ai(t) and 02 (t). Intuitively, we 
have split Eq. (|99|) into two pieces by assigning the mixing term —2A[^^^{t) in equal parts to each 
channels. In the following subsections we show that Eqs. (|100|) - (jl01|) fulfill all the required properties. 

The decay probabilities wi{t) and W2{t), formally defined in Eqs. (f83| and ([84l) . can be easily 
obtained from Eqs. (fT00|) - (fT0T1) as: 

wiit) = Ai(0) + AmAO) - Mt) - A^,,{t) , (102) 
W2{t) - ^2(0) + A„„(0) - A2{t) - A™„(t) . (103) 
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Figure 6: The partial distributions cfg\x) and d^g\x) are plotted and compared to ds{x). The nu- 
merical values of Eq. (IT51) are used. 



In Fig. 7 the functions hi{t) and h2{t) as calculated from Eqs. (|100l) - (|10ip are plotted for the usual 
numerical example of Eq. (|18p and are compared to h(t). Both functions vanish for t = and have a 
qualitative similar behavior to the sum h(t), reaching a maximum for short times and then decreasing 
later on. It is also interesting to notice that both functions show an oscillating behavior on top of the 
exponential law. 

4.3.2 Sum 

The sum hi{t) + h2(t) must obviously give h{t): 

h{t) ^ hi{t) + h2{t) . (104) 

In fact, the quantity h(t)dt represents the probability that the particle S decays between t and t + dt, 
which must be the sum of the probability of decaying into the first channel and into the second channel. 
This property, which indeed has been our starting point leading to our postulate, is easily verified, see 
Eq. dMl) and Eqs. pIin| - pUT|) . 

4.3.3 Limits 32 and gi 

It is important to check that hi{t) and h2{t) fulfill the correct limits when one of the coupling constant 
is sent to zero. In the limit (?2 ^ one has, as expected, that 

hi{t)^~A[{t) ,h2{t)^0 . (105) 

Thus, only the first channel survives and h{t) = hi{t). Similarly, in the limit gi — > one has the 
specular situation 

h2{t) = -A'2{t) ,hi{t)=Q . (106) 
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Figure 7: The partial decay probability densities hi{t) and h^it) are plotted and compared to h{t). 
The numerical values of Eq. (ITS]) are used. 



4.3.4 Integrals and partial decay rates 

The functions hi{t) and h2{t) must fulfill the following requirements: 



«7i(oo)= / dt/ii(t)= / dx4^^(a;) , (107) 



W2(oo) = / dth2{t)= dxd^s\x) . (108) 
Jo Jo 

In fact, u'i(oo) = dthi{t) represents the probability that the particle decays in the first channel 

between (0, oo), a quantity which is also given by dxd!g \x). (This is a typical result of the study 
of line shapes, e.g. Ref. [25l and refs. therein.) 
In order to prove Eq. (|107p we calculate 

) poo 

dthlit) - / dt {~A[{t) ~ A'^,,{t)) = Alio) + A™„(0). (109) 
Jo 

Taking into account that 

^i(O) = (^1°° dxd'^s'\^)) . (110) 

Armm=(yj dxd^s\^)^(^j dx4')(x)^ , (111) 

we obtain: 

/ poo \ / poo poo \ poo 

Ai(0)+A„„(0)= / dx4^(a;)) / dxd^s\x)+ dxd^s\x)] = dxd^s\x) . (112) 

\J —QO / \J — OQ J —OO / J —OO 

In this way Eq. (|107p is proven. The validity of Eq. (|108p can be easily proven following the same 
steps. 
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4.3.5 The Breit-Wigner limit 

Also in the two-channel case it is important to study the Breit-Wigner limit. The (non-relativistic) 
BW distributions d^s-Bwi^) ^^'^ ^s-bw(^) ^'i^- dSHl-dini) are obtained just as in Sec. 2.3: 



S^\x) ^ 4'Rw(a^) - TT- \ . (113) 

s V ; s-BwV ; 27r (a; _ m) V r2/4 ^ ' 



d^^\x) ^ 4'bw(2^) = ^ \ ■ (114) 

Obviously, ds-B\N{x) = 4_bw('^) '^s-bw('^)- "^^^ amplitudes ai{t) and 02 (t) read in the BW limit 

/OO p 
dx4^(x)e-"* = ^e-^*/2e**" . (116) 
-00 ^ 

The quantity Ai, and Amix take the form 



Ai{t) ^ Ai,Bw(i) - |ai3w(i)|' = ( e"" , (117) 



vl2(i)^A2,Bw(i) = |a2,Bw(0r= ( ^1 ^ ' (11^) 
. /.^ . ai,Bw(*)"2,Bw(0 + «i,Bw(0a2,Bw(0 rir2 _rt 

Out of these equations we can easily evaluate the functions hi and /i2 in the BW limit: 

--^^ = Tie"" , (120) 



i2,BWl 







r 


^ r 


^ 2 






r 



= Tae-" . (121) 



The decay probabilities wi (t) and u'2 (t) take the form: 



i«i(i)^^«i,Bw(i) = Y-(l ' (122) 
i«2(i)^w'2,Bw(i)-^(l-e-") . (123) 

In particular, when calculating the decay fractions in a given channel, we obtain the expected results 

f°° Ti 
wi,Bw(oo) = / dtft,i_Bw(0 = ' (124) 

/■°° r2 

W2,Bw(oo) = / dtft,2,Bw(0 = -fT • (125) 

In Fig. 8 the functions hi{t) and h2{t) of Eqs. (|100p - (jlOip are plotted separately and compared to 
the corresponding exponential functions /ii.Bw(i) smd /i2,Bw(0- It is visible that also in each channel 
the deviations from the BW limit are not negligible for small times, but become smaller and smaller 
with increasing t. It is also interesting to notice that the two functions hi{t) and h2{t) oscillate around 
the corresponding BW limit in a different and peculiar way. 
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t [2m2]"'' t [Im^] 

Figure 8: Left: the decay probability density into the first channel hi{t) is plotted together with the 
BW limit /ii.Bw(i) = Tie"^*. Right: the same for hi{t) and ft.2.Bw(i)- In both cases the numerical 
values of Eq. have been used. 



The results obtained in the BW limit are expected: the usual exponential law is recovered and the 
different channels are proportional to the corresponding partial decay width T^. In the BW limit the 
ratio R{t) defined in Eq. (|57)) becomes a constant: 

^(^) = r4S ^ RBw{t) = ^ = const . (126) 

''■2{t) 1-2 

This result does not hold true in the general case, in which R{t) varies with t. In Fig. 9 the function 
R{t) and the constant BW limit i?Bw(i) = ri/r2 are shown: the deviation of the exact result R{t) 
from the BW limit is sizable. There are time intervals in which the production of (pi(pi pairs is 
enhanced with respect to ip2^2 pairs, and there are time intervals in which the opposite is true. This 
is an extremely interesting behavior which shall be studied in detail in the future. The reason for the 
large departure from the constant BW result can be also traced back to the oscillating behavior of 
both functions hi{t) and h2{t). The frequencies of the oscillations are different, and thus the ratio R{t) 
results as a complicated outcome. Moreover, the deviations are still large also for large times. This is 
due to the fact that for t 3> tbw the quantities hi(t) and h2{t) are both small numbers: even small 
variations from the exponential limit cause large deviations of the ratio R{t) from the constant value 
ri/r2. The full study of the large time behavior of the function R{t) represents also an interesting 
outlook for future studies. 

Finally, in Fig. 10 we plot the ratios 

Both functions tend to the constant value Ti/T in the BW limit. One can see that hi{t)/h{t) oscillate 
and that Wi{t)/w{t) tend to a constant limit which arc not coincident with the naive BW expectations, 
thus showing that deviations persist also for large times. 

4.3.6 Invariance under time-reversal 

The underlying theory described by the Lagrangian of Eq. ^ is invariant under time reversal; there- 
fore, the physical quantities under study should also reflect this basic symmetry. The amplitude a{t) 
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R(t) 




Figure 9: The ratio R{t) ~ hi{t)/h2{t) (solid line) is plotted and compared to the constant BW limit 
-RBw(i) — ri/r2 (dashed line) Large deviations, which persist also at large times, take place. As 
usual, the numerical values ([T8| have been used. 




Figure 10: Left: hi{t)/h{t) and wi{t)/w{t) have been plotted and compared to the constant BW limit 
Fi/r. Right: h2{t)/h{t) and W2{t)/w{t) have been plotted and compared to the BW limit F2/F. The 
numerical values of Eq. (|18p have been used. 
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defined in Eq. (fT9|) fulfills the equation a{—t) — a*{t). As a consequence, the survival probability 
p{t) = a*(t)a(t) is invariant under time reversal: p{—t) ~p{t). 

Also the decay probability w(t) in Eq. (j24p is invariant under time reversal for the very same 
reason: w{—t) — w{t). Notice that the decay probability density h{t) is odd under time reversal: 
h{—t) = —h{t). In this way the invariance of the physical quantity h{t)dt, and therefore w(t), is 
assured. 

The invariance under time reversal is satisfied also by the decay probabilities wi(t) and W2(t), 
as the formal expressions in Eqs. ((55| and ([M|) show: ■wi{—t) = wi{t), W2{^t) — W2{t). Due to 
the mathematical property ai.2(— = a* 2(^) ^-Iso the expressions for wi{t) and W2{t) in Eqs. (I102p 
and (|103p fulfill the requirement of invariance under time reversal. Just as for the function h{t), the 
decay probability densities hi{t) and h2{t) must be odd under time reversal: hi{—t) = —hi{t) and 
h2{-t) = -h2{t). This is easily verified for Eqs. (fTOO|) and (fTOTj) . 

4.3.7 On the uniqueness of the postulated solution 

Eqs. (|100p - (|101l) fulfill all the required conditions to represent the decay probability densities in the 
first and the second channel. However, it is important to study the uniqueness of presented solutions. 

A different combination hi{t) = -A'^{t) - aA'^i^{t), h2{t) = -A'^it) - (3A'^i^{t) with a + /3 = 2 
and a, (3 ^ 1 would obviously satisfy the constraint given by the sum (|104l) (Sec. 4.3.2) and the limits 
91,92 — >■ (Sec. 4.3.3), but would not fulfill the constraints given by the integrals (|107p - (|108l) and 
would not reproduce the correct BW limit. Therefore, only the case a = (3 = 1 is allowed. 

More in general, one can try to find alternative solutions by writing down the functions 

hi{t) = hi{t) + 5{t) (127) 
hit) ^ h2{t) ~ S{t) (128) 

where the 'hypothetical' function 5{t) = 5{gi,nii, g2,rn2,t) has been introduced. While the sum 
constraint of Eq. (|104p is fulfilled by construction, the function S{gi,mi, g2,m2,t) must vanish for 
gi — > and for g2 ^ (Sec. 4.3.3) and the integral over time must vanish, 5{t)dt ~ 0, in order 
that Eqs. (|107p - (|108p still hold. Then, 5{t) must also vanish in the BW limit to be in agreement 
with Sec. 4.3.5. Moreover, being hi{t) and h2{t) positive definite, \5{t)\ < h2{t) when S{t) > and 
\6{t)\ < hi{t) when 5{t) < 0. All these properties are automatically fulfilled by the trivial solution 
S{t) = 0, but the question is if other possibilities exist. 

A further property that S must fulfill is obtained by studying the exchange of the two particles. 
This operation amounts to exchange gi < — > 52 and mi — > m2. While in this case fti /12, it must 
hold that: 

S{gi, mi, g2, 1712,1) = S{g2,m2, gi,mi,t) . (129) 
If we consider, for instance, the function 

6it) = do^lTn[aiit)a*2it)] (130) 
at 

where the constant Sq is small enough to guarantee the semi-positive definiteness of hi{t) and h2{t), 
it is easy to prove that Eq. (|130p would actually fulfill all the mentioned properties from Sec. 4.3.2 
to Sec. 4.3.5 together with the exchange property of Eq. (|129p . However, it violates time-reversal, 
because it is even under switch of t: S{—t) = 6{t) (while, as shown in Sec. 4.3.6, the decay probability 
densities must be odd under time reversal). For this reason the only possibility is 60 — 0, i.e. the 
function d{t) vanishes. 

This discussion can be generalized. Writing S{t) = ^A(t) and making the additional assumption 
that the function A{t) can be expressed as A{t) — "Ylii j Cija,i{t)a*(t) where the coefficients Cy are pure 
numbers, it is not possible to build 8(1) such that it is odd under time reversal and at the same time 
fulfills Eq. (jl29l) . This chain of arguments represent a partial proof of the validity of Eqs. (I100p - (|10ip . 
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Summarizing, all these arguments enforce our point of view that the Eqs. ()100|) - (|101|) . and therefore 
Eqs. (|102l) - (|103p . are the correct answer to the problem. However, the final answer can only be given 
if one can mathematically prove the validity of the equations (|100p - (|10ip under general assumptions. 
Presently, we did not achieve this point, which remains an important task for future studies. However, 
by using the formalism of the Lee Hamiltonians it is possible to find the exact expressions for hi(t) 
and Wi{t), which allow for a numerical verification of Eqs. ()100p - (|10ip . see Sec. 5. 

4.3.8 Generalizations 

The present discussion can be easily generalized to the case of a n-channel decay: this case is discussed 
in Appendix B. Moreover, the formalism is not restricted to scalar field theories only but can be easily 
generalized to each QFT Lagrangian. 

4.3.9 Deviations from the exponential decay law and measurements 

In this work we focus our attention to the undisturbed decay probabilities of the system: after the 
creation at t — the unstable system \S) evolves according to the Hamiltonian of the system which 
control its decay. We do not attempt here to include the measurement into the mathematical treatment. 
This is surely an interesting subject which is left as an outlook. 

This problem has been first investigated in Refs. [Tl|29l|30]. In particular, in Ref. [29] it has been 
pointed out that random measurements occurring with a certain mean frequency v make the decay still 
appear to be an exponential, even if the frequency is large enough to probe the non-exponential decay 
time. This is indeed a generalization of the simple discussion of Sec. 2. 5, where the quantum-Zeno 
effect has been presented as an outcome of pulsed measurements. It is assumed that the measurement 
lets the state collapse back into the unstable state \S), which also implies a reset of the clock to the 
initial value of i = 0. In this framework, the deviations from the exponential decay law studied in this 
work would manifest in an exponential decay laws with a ratio which can be different from ri/r2 in 
dependence of the (mean) frequency of the observations. 

On the other hand, future studies in this direction should also take into account the theory of 
measurement, see Ref. [3TI and refs. therein, in which it is shown that the influence of the measurement 
apparatus and its response time can have a non-negligible effect. As studied also in Refs. [32l|35, the 
outcome may also depend on the details of the experimental apparatus. A detailed investigation of the 
'collapse' beyond the simple projection postulate is also an interesting direction of study in connection 
with non-exponential decay laws. Another interesting line of research is to investigate in more details 
the problem of sequential decay processes analyzed in Ref. [30] by using the approach described in the 
present manuscript. 

5 Two-channel case: QM and link to QFT 

5.1 Exact determination of the decay probabilities Wi{t) and W2{t) in QM 
by using the nonrelativistic LH 

We now study the two-channel decay problem in the QM framework by using the Lee Hamiltonian H 
introduced in Eq. (|49|). 

Due to the formal analogy of the QM-LH and the QFT approaches, we can repeat here all the 
steps done in Sec. 4.2 and 4.3 leading to the postulated solutions for Wi{t) and hi{t). They read: 



W^{t) = A,(0) + A,mx{0) - A^{t) - Armx{t) , h,{t) 

A,{t) ^\a^{t)\^ , A„,^At) = Rc [ai{t)a*2{t)] , 



(131) 
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where 



dxd^ (a;)e 



.gflm[S,(x)] 



(x - Mo + Re[n(2;)])2 + W [U{x)] 



(132) 



(See Sec. 3.1 for the definitfons of the functfons Si (a;) and n(.T) m the nonrelativistic QM-LH case). 

In the present QM-LH treatment we can calculate exactly the functions Wi{t). To this end we first 
write the formal expression for the decay probability 'Wi{t) in the i-th channel: 



^{t)= / /i,(M)du = ^|(i,k|e-*-^*|S')f ,i = l,2. 



We thus need to evaluate the matrix clement 



(z,k| 



-iHt I 



s) 



Proceeding as in Sec. 3.1 we obtain the following expression: 



(z,k| 



-iHt\ 



V J- 



dE 



Gs{E) 



E — Wj(k) + is 



(133) 



(134) 



(135) 



Introducing the dispersion relation of Gs{E), see Eq. 
we arrive at the following equation: 



67|) . and performing the subsequent integral, 



(*,k| 



iHtl 



V J- 



dxds{x) 



Wi(k) 



Using now Eq. (jl33p we obtain the final, exact expression for Wi{t) (in the continuos limit) as: 



(M 



dxds{x)- 



-iuJi {k)i 



, /if ^'='(t) = -- 



dt 



(136) 



(137) 



We have now compared the expressions in Eq. (|137p and Eqs. (|13ip - ()132p for different parametriza- 
tions of the functions a;i(k) and /i(k). In particular, we have used: 

(i) A generalized version of the two-pole model developed in Ref. [1], in which two open channels 
are included. The numerical values have been varied in order to check dependence on them. 

(ii) Modifications of the two-pole model including energy thresholds. 

(iii) Different choices for the loop functions /i(k). 

In all these cases we have verified that the two expressions ()137|) and (|131|) - (jl32|) are numerically 
indistinguishable. All the results obtained in the present QM-LH approach lead us to conclude that the 
postulated expressions in Eqs. (|13ip - (|132l) (which have the same form as the QFT-like expressions in 
Sec. 4.3) represent -at least- a very good approximation of the exact results. A proof of the complete 
equivalence (or eventually of the -albeit small- deviations) of these expressions is left for the future. 
In particular, one should check the validity of the Ansatz in Eqs. ()131|) - (|132|) under more complicated 
circumstances, involving a numerical check over a large range of parameters and with different forms 
of the functions /i(k). 

In the end, it is also important to stress that the qualitative behavior of the functions hi{t) and R{t) 
as evaluated in the present QM-LH approach is very similar to the plots presented in this manuscript: 
hi{t) vanish for t — and then tend to the BW limit for large times. Typical oscillations are present 
and the ratio R{t) is varying abruptly as in Fig. 9. 
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5.2 Generalization to the relativistic case and link to QFT 

We now study the two-channel problem using the formalism developed in Sec. 3.2, in which a relativistic 
generalization of the QM-LH approach was developed. 

The solutions postulated in Sec. 4.2 and 4.3 takes the following form in the present case: 



Mt) = \a,{t)f 



mix 

A 



(0) - A,{t) - Arm At) , h,{t) = -w[{t) , 
,{t)^Re[ai{t)a;{t)] , 



where 



ai{t) 



- M§ + Re[n(a;)])2 + W [U{x)] 



(138) 



(139) 



(See Sec. 3.2 for the definitions of the functions I],;(x) and n(x) in the relativistic LH treatment). 
Note that, when using the identifications of Eqs. ([7^ - (|M|) . the previous equations coincide with those 
in Sec. 4.3. 

By following the same steps of Sec. 5.1 (that is, taking into account the relativistic form of the 
propagators) we arrive at the following 'exact' expressions for Wi{t) for the relativistic LH approach: 



Also in this case we have numerically compared the conjectured solutions of Eqs. (|138p - (jl39p with the 
exact result in Eq. (jl40|) and found them equal. 

In particular, when using the identifications of Eqs. (|79 |) - ([80l) which lead to a mathematical equiva- 
lence of the present LH approach with the QFT model described in Sees. 2 and 4, we have recalculated 
all the numerical results presented in Figs. 2-10 with the help of the exact expressions in Eq. (|140l) : as 
expected, no deviation was found. This fact reinforces our point of view that the expressions represent 
the exact result (or a very accurate approximation of it). 



dxds{x)- 



di 



(140) 



6 Conclusions 

In this work we have studied the decay law(s) as function of time of an unstable particle S, which can 
decay into two (or more) decay channels. We have done this investigation in two different theoretical 
contexts: using a relativistic quantum field theoretical Lagrangian and using a quantum mechanical 
approach expressed in terms of a Lee Hamiltonian. 

In the first part of the manuscript we have reviewed general properties of the survival probability 
p{t) and the decay probability density h{t) = —p'{t) (Fig. 1-5, Sec. 2 and Sec. 3) for both the QFT 
and the QM cases: the non-exponential behavior of the survival probability p{t) for short times, the 
Breit-Wigner limit leading to the usual exponential decay law, and the quantum Zeno effect have been 
described. Formally, the QFT and the QM-LH approaches are very similar, the main difference being 
the presence of relativistic propagator in the former case and nonrelativistic propagators in the latter 
one. To overcome this difference we have also developed a simple relativistic generalization of the 
QM-LH approach in which the same equations of the QFT case (can) arise. 

In the second part of the manuscript (Sees. 4 and 5) we have concentrated our attention on the 
peculiarities of the decay into each one of the two (or more) available open channels. In the two- 
channel case the mathematical expressions for the decay probability densities into each channel, hi(t) 
and h2{t), and of their corresponding integrated quantities wi{t) = du/ii(u) and W2{t) = Jq du/i2(u) 
have been discussed. Namely, hi{t)dt represents the probability that the unstable particle S decays 
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into the i-th channel between t and t + dt, while Wi{t) represent the probability that the particle decays 
into the i-th channel between and t. Obviously, wi{t) + W2{t) — 1 — p{t). 

In particular, in QFT (Sec. 4) we did not evaluate hi{t) and Wi{t) exactly, but we have postulated 
a solution (Eqs. (|100l) - (|10ip in Sec. 4.3, see also Figs. 7-8). We have shown that the proposed 
expressions for hi{t) and Wi{t) fulfill all the required properties, but we could prove the mathematical 
exactness of the results only under a restricted set of hypotheses, see Sec. 4.3.7 and below. We have 
also shown that the conjectured solutions show strong deviations from the Breit-Wigner limit, in which 
the decay law is exponential. An important quantity that we have studied is the ratio of probability 
densities R{t) — hi{t)/h2{t). While this quantity is constant in the Breit-Wigner limit and equals the 
ratio of the tree-level decay widths {R{t) — ri/r2), it turns out that the full expression is in general 
not constant and may deviate sizably from the constant limit, see Fig. 9. The functions hi{t)/h{t) and 
Wi{t)/w{t) show also the peculiar characteristics of the two-channel problem and have been plotted in 
Fig. 10. 

In view of the fact that the conjectured solutions in QFT could not be analytically demonstrated 
in the general case, we have studied the same problem in the formally very similar (up to the presence 
of nonrelativistic propagators) approach of a nonrelativistic QM-LH approach. In this case, besides 
the postulated solutions, it is also possible to derive the exact form for the functions hf^°-'^*{t) and 
^exact(^j.y numerical analysis was performed in order to show that the proposed solutions hi{t) 
and the exact ones hf^°^'^^{t) coincide, thus proving that the proposed solutions represent -at least- a 
very good approximation of the exact results. When studying the relativistic generalization of the Lee 
Hamiltonian approach, in which the equations can be chosen to coincide with the QFT study, we have 
once more numerically verified that the expressions postulated in the QFT context are valid: using 
the very same parameters of Eq. we have checked the correctness of Figs. 2-10. 

The are various possible outlooks for the future, which can be divided into mathematical develop- 
ments and physical applications. 

On the mathematical side, one can perform the following studies: (i) Mathematical proof of the 
correctness of the postulated expressions for hi{t) and 'Wi{t) presented in Sec. 4.3 in the QFT case. If 
the proposed solutions are not exact, it would be necessary to investigate why they offer such a good 
numerical approximation to the problem. In this context, one needs to check the numerical validity of 
the postulated Ansatz under more general conditions, (ii) Inclusion of higher orders in the calculation 
of the self-energy diagrams of the QFT case, thus going beyond Fig. l.b. (iii) Extension to (strictly) 
renormalizable Lagrangians. 

On a phenomenological level various investigations are possible: (i) Study of physical processes 
in the framework of the Standard Model, where decays of unstable particles into many channels have 
been observed and precisely measured. To this end the mathematical development of the mathematical 
outlook (iii) is necessary, (ii) Applications to weak decays (such as electronic capture decays) of 
atoms. This analysis could be important to study in detail the properties of oscillations such as 
the ones experimentally found in Ref. [26], see also Ref. [11] for a related theoretical discussion, 
(iii) Applications to hadron physics: the deviations from the non-exponential law are large in this 
case [T^, although they last very short due to the short- living nature of (most) hadrons. However, 
deviations from the exponential decay law might be interesting in connection to an expanding gas of 
hadrons produced in heavy ion collisions [27] . (iv) A peculiar situation takes place when one channel 
is kinematically not allowed, M < 2m2- When M is just slightly smaller then 2m2, the finite width of 
the unstable state S due to its decay into the first channel allows also the decay into the -nominally 
forbidden- second one. This situation is physically realized in the case of the scalar mesons ao(980) 
and /o(980), see details in Refs. [TH [T9||25l [28] and refs. therein. The study of the temporal evolution 
in such a case represents an interesting outlook. 

As a last remark, it is conceivable that the decay into different channels can be experimentally 
studied with the procedure of Ref. [8] ; asymmetric barriers should be used, in order to allow for different 
tunneling probabilities is different directions. Then, it would be possible to test in a experimentally 
controlled way quantities such as the ratio R{t) = hi{t)/h2{t) described in this work. More in general, 
in connection with the experiments, one should also go beyond the undisturbed time-evolution of the 
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system and include the measuring procedure as a part of the physical description of the decay. 

Acknowledgments: The author thanks Giuseppe Pagliara for cooperation and many valuable 
discussions on the subject and Thomas Wolkanowski for useful remarks. 



A The loop function 

In order to regularize the self-energy diagram of Fig. l.b we introduce at each vertex the vertex- function 
(t>{q). At the level of the Lagrangian formulation this can be achieved by rendering the Lagrangian 
nonlocal, see details in Refs. (THl [M] and refs. therein. The integral in Eq. pU)) takes the form 



m-' + ie 



(141) 



A general property of E(a;, m^) follows from the optical theorem: 



(V25)2lm[I](x,m2)] = xT'-\x,m,g) U{q = (0,q)) 



where 



(142) 



Due to spatial isotropy, the function (t){q = (0, q)) must depend on q^. 

For a generic (/'(q), the tree- level decay width is indeed not given by Eq. ([2]) but takes the following 
modified form: 

T'-\x, m, g) ^ T'-\x, m, g) [^[q = (0, q))] ' . (143) 

When assuming that 0(g) — 0(q) (i.e. when we drop the dependence on g^), we can perform, in 
the rest frame of the S particle [p ~ {x, 0)), the integral over q'^ in Eq. (|14ip . obtaining: 



I](a;2,m2) = 



(27r)3 ^q2 + rn'^ (4(q2 „ 2-2 „ 



In this work we make the following simple choice for numerical evaluations: 

m - 0(q) - ^(A' - q') , 



(144) 



(145) 



i.e. we work with a sharp cutoff A. (In general, the use of smooth functions affects the results only 
slightly, thus showing a weak dependence on the adopted vertex function. An explicit calculation has 
been presented in Ref. [15]). In this case, by performing also the spatial integration we obtain the 
explicit result: 



^{x, m^) = 



— V 4m^ — , 
8n^x 



■ arctan 



Ax 



VA^ + m-^^/Am'^ — x'^ 



1 

8^ 



log 



A -I- VA2 + m2 



(146) 



Following comments are in order: 

(i) The imaginary part of the self-energy amplitude Im[E(a;^, m^)] is zero for < x < 2m and 
nonzero starting at threshold. 

(ii) The real part {\/2g)'^Ke[T,{x^ , m^)] is nonzero below and above threshold and depends explicitly 
on the cutoff A. 
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(iii) For the case ^(q) = 6{h? — q^), the optical theorem takes the form 



{V2gflwL[Y.{x^ ,w?)] = xT^-\x, m, g)e + m? - |) , (147) 

Besides the additional function, it does not depend on the cutoff A. Note that Eq. (fT3|) in the text 
is strictly speaking valid only for x < 2y^A^"+mf. This condition is numerically met in Fig. 2 and 
Fig. 6. 

(iv) The present choice of 4>{q) breaks Lorentz invariance. This is here not a problem, because we 
always work in the rest frame of S. It is also not difficult to generalize 0(q) in order that it is Lorentz 
invariant and delivers the same results of this paper when the rest frame of S is considered. 

As a last step we turn to the issue of the bare and renormalized masses: in this work we started 
with the bare mass Mq entering in the Lagrangian of Eq. ([1]) and then we derived the renormalized 
mass M in Eq. (jl6p . arising upon the inclusion and the resummation of the self-energy diagram of 
Fig. l.b. Alternatively, one could have added a counterterm in Eq. ([1]): 

£^ C-^CS"^ with C = R{Ma). (148) 
In this way the mass equation takes the form 

- Ml -C + R{M) = 0^ M = Mo . 
Clearly, all the physical results of this work would be unaffected by this alternative procedure. 



B The n-channel case 

In the n-channel case the interaction Lagrangian takes the form 

n 

Ant =^5.5^', (149) 

1=1 

where the particle ipi has a mass m^. The spectral function (or mass distribution) ds{x) can be 
decomposed as ds{x) — X]"=i d-s\x) with 

rf«(x) = ^ hm 2gflm[S(pW)]+e ^ ±iV2g.mx\mj) . (150) 

^ e^o (a.2_A^2^Ren(x2))' + (Imn(x2)+e)' ^ ^ l-A^ > ^ ^ ^ ' 

We define [t) as the Fourier transform of cfg {x) : 

/oo 
dx4^(a;)e-"* (151) 
-oo 

Out of ai(t) we define 

„ " a,(t)a*(t) + a*(t)ai(t) 

The probability that the state \S) decays in the time interval {t,t + dt) into the i-th channel reads 
hi{t)dt, where 

h,{t)^-A[{t)-A'„,,,^,{t). (153) 
It is possible to define the n{n~ l)/2 ratio functions 

= TtI ^ith z < J, I - J = 1, ...n. (154) 
hj{t) 
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In the BW limit one has Rij (t) — Ti /Tj , where Ti is the tree-level decay width in the i-th channel. In 
general, we expect large deviations from the BW limit. Indeed, many unstable particles of the Standard 
Model decay in more than two channels: the analysis of the temporal behavior can be performed 
through the formulae derived here. Explicit calculations represent an outlook for the future. 
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